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Abstract 

We study the associated production of a charged Higgs and a W boson in high-energy pp cohi- 
sions at the Large Hadron CoUider. This is an interesting process for charged Higgs discovery, or 
exclusion, since the production cross section could depend strongly on the model, offering poten- 
tial discriminating power between super symmetric extensions of the Standard Model with minimal 
or extended Higgs sectors. We compute the cross section for this process in the next-to-minimal 
supersymmetric Standard Model (NMSSM), at the tree level for quark-quark scattering and at 
one-loop level for gluon-gluon scattering. The most important corrections beyond leading order 
are taken into account using an improved Born approximation. We find that the pp — t- H^W^ 
cross section can be resonantly enhanced by up to an order of magnitude over its MSSM value 
(both for ^/s = 7 TeV and 14 TeV) through the contributions of heavy, neutral, singlet-dominated 
Higgs bosons appearing in the s-channel. Since such Higgs mass configurations are normally not 
possible in the MSSM, the observation of associated H^W^ production at the LHC could provide 
a striking, although indirect, signature of a more complicated Higgs sector. 
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I. INTRODUCTION 



There are good phenomenological and theoretical reasons to consider supersymmetric 
(SUSY) extensions of the Standard Model, including amelioration of the hierarchy problem 
or the fine-tuning of the Higgs mass, unification of the gauge couplings, and the existence 
of a dark matter WIMP candidate. The most common SUSY model is the minimal super- 
symmetric extension of the Standard Model (MSSM) [H, 0], which exhibits all of the above 
benefits. These benefits come at the price of introducing a mass term (called the //-term) 
for the two Higgs doublets that are required in the MSSM. This leads to another fine-tuning 
problem, namely why the mass parameter /i of the Higgs doublets should be at the elec- 
troweak mass scale, as is required by phenomenology. This is known as the /i-problem. 
Furthermore, the mass of the lightest Higgs boson cannot be less than ~ 115 GeV, while 
in the MSSM at the tree-level it has an upper limit equal to the Z boson mass. To fulfill 
the experimental constraints, large corrections to the Higgs mass from top- and stop-loops 
are needed, which require a rather large stop mass leading to an additional amount of fine- 
tuning in the model. If we take fine-tuning arguments seriously, we may therefore consider 
whether there are alternatives to the MSSM (see e.g. js], 0]). 

One such alternative is the next-to- minimal supersymmetric Standard Model (NMSSM), 
which has recently been reviewed in Refs. js], [6|. In the NMSSM, there is in addition to the 
two Higgs doublets of the MSSM a singlet Higgs field, which is the scalar component of a 
chiral singlet superfield added to the MSSM superpotential. The reason for introducing this 
additional scalar is that the /i-term is now dynamically generated, so that the fine-tuned 
parameter /i is no longer needed. This scalar mixes with the other scalars from the two 
doublets, leading to a Higgs sector with seven Higgs bosons, compared to the five present 
in the MSSM. The fermion component of the singlet superfield, the singlino, additionally 
mixes with the neutralinos, providing interesting possibilities for dark matter that can be, 
e.g., singlino or singlino-Higgsino-dominated jsj, or very light 0, Isl. 

There are two extra neutral Higgs bosons in the NMSSM compared to the MSSM; one 
CP-even and one CP-odd. The tree-level mass relations for the Higgs bosons are then also 
modified, and it is possible for one or more Higgs bosons to be quite light. In particular the 
lightest CP-odd Higgs, Ai, can be significantly lighter than in the MSSM — viable scenarios 
with < 2mfe exist joj. The charged Higgs boson can also be rather light [lo|, 11|, albeit 
not as light as the A\. 

Because of the different Higgs phenomenology due to the modified mass relations and the 
additional particles, the parameters of the model are not as constrained as in the MSSM, 
and new decay channels and production mechanisms may become important at the LHC. 
For example, the charged Higgs may decay as — W~^Ai, and the Ai in turn may 
decay dominantly as Ai bb, t^t^ . Such differences are important to take into account in 
searches at LHC, so that no possibilities are missed. 

In this paper, our focus is on production of the charged Higgs boson in association 
with a boson. This is not the main production channel usually considered for and, 
to the best of our knowledge, it has not been studied in the context of the NMSSM. It 
was however pointed out in Ref. T3] that a related process, the associated central exclusive 
production of H^W^, may be useful in the NMSSM. 

The usual production mechanisms for are top quark decays t — )■ H~^b for light 
charged Higgs (where the charged Higgs boson mass m^i < rrit), and production with 
a top, bg — )■ H~t or gg — H~tb, for heavy charged Higgs {mH± > mt). The associated pro- 
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duction mechanism may, however, be important to secure additional information about the 
Higgs sector if the charged Higgs is first observed in one of the above-mentioned processes. 
An advantage of the H'^W^ process is that a leptonically decaying W may be used as an ex- 
perimental requirement. The cross section for associated production is as we shall see rather 
mo del- dependent, and the observation of this process may therefore provide constraints on 
the model parameters. For example, in the NMSSM, smaller tan/3 is allowed, leading to 
a possible enhancement of the cross section. A potentially more important difference be- 
tween the MSSM results and the NMSSM comes from the resonant s-channel exchanges 
of additional singlet-dominated Higgs bosons. Due to the very restricted mass relations of 
the MSSM Higgs sector, these contributions cannot be resonant in the MSSM, while in the 
NMSSM they can. The resonant enhancement of the parton-level cross section will also en- 
hance the hadron-level cross section in some range of charged Higgs boson masses and this 
enhancement could be potentially visible at the LHC It could be important for discerning 
differences between MSSM and NMSSM, and for setting limits on the parameter space. 

The H^W^ production channel was first considered in [l3|, where the cross sections were 
calculated in the approximation that = 0. This study did not include any contribution 
from squark loops. The cross sections were later calculated in full generality for two-Higgs- 
doublet models (2HDM) and the MSSM in 14-17|, where in particular [l6|, 17| included the 
squark loop contributions. In [isj, it was further shown that there may be a substantial 
enhancement of the cross section compared to the MSSM in a general 2HDM. Our paper 
extends these studies to the NMSSM, and we are going to investigate in particular the 
differences between the MSSM and the NMSSM for this process. 

Next-to-leading order (NLO) corrections to the bb — )• H'^W^ subprocess, including QCD, 

SUSY QCD, and electroweak contributions, are known for the MSSM [iSlliL 

There have 



also been phenomenological studies of this process for the LHC at 14 TeV 0, • While 

we do discuss the LHC aspects of the H^W^ process in the NMSSM below, the inclusion of 
NLO corrections and a detailed study of LHC signatures are beyond the scope of this paper 
and will be left for future studies. 

Finally, we restrict ourselves to the case of charged Higgs bosons heavier than the top 
quark, which means that the t — )■ bH~^ production mechanism is not effective. The reason is 
twofold: first, this decay is the same in the MSSM and the NMSSM, but we are interested 
in differences between the two models. Second, the experimental constraints on this decay 



channel are already quite strict (29|-l3ll, ruhng out BR(t bH+) > 5-10% (in the MSSM). 



II. THE NMSSM HIGGS SECTOR 

The NMSSM is defined by removing the yU-term from the MSSM superpotential and 
adding a singlet chiral superfield S, which only couples to the Higgs doublets. Assuming 
scale invariance, the general form of the superpotential is 



W^vfMSSM = W^MSSM + XSHuHd + 2 



(2.1) 



In the above superpotential, the (unchanged) Yukawa terms are contained in VFmssm- The 
scalar potential of the NMSSM is obtained from the F- and D-terms plus the soft SUSY- 
breaking terms for the Higgs sector. 



soft 



\Hu 
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XAxSHuHd + -KA^S' 



h.c. 



(2.2) 



3 



where the dimensionless couphngs A and k, the soft SUSY-breaking parameters Ax and 
Ai^ with dimension of mass, and the singlet mass rris are new parameters compared to 
the MSSM. As usual, rris is fixed by the minimization of the potential. Requiring that S 
acquires a vacuum expectation value (vev), s = {S), yields an additional new parameter of 
the model, and gives rise to an effective //-parameter fi^s = As. Together with the ratio of 
vevs of the two Higgs doublets, tan/3 = v^/vd, where v"^ + v"^ = v"^ = (174 GeV)^, we have 
six free parameters of the Higgs sector of the NMSSM: A, k, Ax, A^, s, and tan /3. 

As S* is a complex field, there will be two additional physical Higgs bosons in the NMSSM 
compared to the MSSM. For a CP-conserving theory (as is assumed here) we have three 
CP-even neutral states Hi,H2, H3 and two CP-odd neutral states Ai and A2, where we take 
the states to be ordered in mass with Hi and Ai the lightest states. 

The mass of the charged Higgs boson is at tree-level 

= ^^(v4A + Ks) + m^-AV, (2.3) 

which can be compared with the MSSM expression 7Ti^± = m^ + m^. It therefore simplifies 
our expressions to define an effective "doublet mass" in the NMSSM as 

< = + ^S) = ^{Ax + KS) = ^{Ax + ^/.eff), (2.4) 

sm/p spCjs spcp A 

where we defined sp = sin/3 and C/3 = cos/3. Thus, for fixed fics, the squared doublet mass 
and the charged Higgs mass both depend linearly on Ax and on k/X. The parameter Ax 
may be everywhere swapped for tua as a parameter of the Higgs sector. The mass relation 
(12. 3 P can then be written as 

m^± = m\ + — A^f ^. (2.5) 

This exhibits one important difference between the NMSSM and the MSSM: in the MSSM, 
the charged Higgs and the CP-odd Higgs are almost degenerate in mass as soon as they 
are heavier than ~ 200 GeV. In the NMSSM there is no such strong correlation; partly 
because of the additional contribution —X^v"^ which lowers ni?^±, but mainly because there 
is usually no physical state with mass niA- The two CP-odd states in the NMSSM arise as 
mixtures of the CP-odd state of the MSSM, A = Im(ifuC^ -|- H^Sis), and of the imaginary 
part of the scalar S. The effective mass thus only corresponds to a physical mass if the 
mixing between the two pseudoscalar bosons vanishes. This mixing is obtained from the 
mass matrix for Ai, A2, 

(m\ — {m\spCp — 3Aks^) \ 

' . (2.6) 

- {m\spCii — 3Aks^) — {m\si3C[s + 3Aks^) — SkA^s J 

In terms of the weak basis eigenstates AJ^^^ = [ImHd, ImHu, ImS*), the physical CP-odd 
eigenstates A'^'^^^ = {Ai, A2) (ordered in increasing mass) are given by A^'^^^ = PijAj*^^^ 
with the 2x3 mixing matrix Pij, or explicitly 

Al\ ^ fPu P12 Pl3 

A2J \-p21 P22 -P23 
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It has been argued, for example by Dermisek and Gunion jo], |32-34|, that Ai may be much 
hghter than the other Higgs bosons, and can even be as hght as a few GeV, and still be 
allowed by EWPT and collider constraints. 

The masses of the CP-even Higgs states require a three-dimensional mixing matrix Stj 
rotating the weak basis HJ'^^^ = (Reif^, Reif^j, ReS") to the physical one, 

rjmass o Tjweak 
Hi = bij^j 

such that the physical mass eigenstates H^^^ are ordered in increasing mass. The corre- 
sponding expressions for masses and mixings are omitted here; they may be found in e.g. jsf. 
It can be shown that the tree-level mass of the lightest Higgs boson Hi is no longer limited 
by mz as in the MSSM, but instead by 

rriH^ < m% cos^ 2/3 + A^;^ sin^ 2/3. (2.7) 

To summarize, in the NMSSM (at tree- level), the lightest CP-even Higgs is allowed to be 
somewhat heavier, and the charged Higgs somewhat lighter, than in the MSSM, while the 
lightest CP-even Higgs may be much lighter than in the MSSM. 

In the above discussion we have considered only the tree-level masses, but just as in the 



MSSM there can be considerable corrections to these masses at higher orders |35H42| . To 



take these into account to the best precision available [43|, we use the code NMSSMToOLS 



(version 2.3.5) [4J, |45| for our numerical evaluation of the Higgs mass spectrum. We also 
use this code in the following to calculate the mixing, all coupling strengths, and the Higgs 
decay widths from the given input parameters. 

Beyond leading order, the Higgs spectrum depends on all the parameters listed above, 
as well as on the soft SUSY-breaking parameters of other sectors; the most important 
corrections typically come from stop mixing. The standard way to cope with this situation 



in the MSSM is to consider a benchmark scenario (such as those defined in [46|, |47j) to 
fix the higher order corrections from other SUSY sectors, and then vary independently 
the parameters in the Higgs sector. We shall use the same approach here, extending in a 
straightforward manner the MSSM benchmark scenario to the NMSSM. We therefore use 
/Xeff as an input, which together with the NMSSM coupling A determines the value of the 
singlet vev s. Since we are mainly interested in comparing the NMSSM to the MSSM, 
rather than the MSSM to itself, we use the same benchmark (inspired by maximal mixing) 
throughout this work: 

MsusY = 1 TeV, = A - fi,s cot /3 = V6Msusy, = = At 

fies = 250 GeV, Mi = 100 GeV, Ma = 200 GeV, M3 = 1 TeV 

For the remaining (free) input parameters of the NMSSM Higgs sector, we shall take Ax, 
tan/3 — which are equivalent to the two parameters mH±, tan/3 in the MSSM through 
Eqs. (12. 4p . (12. 5p — and in addition A, k, and A^^, which are specific to NMSSM. For any 
choice of NMSSM parameters, the corresponding MSSM limit can be obtained by taking 
A — 7- 0, K — 0, while keeping the ratio n/X and all dimensionful parameters fixed. 
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FIG. 1. Leading-order diagrams for the quark-initiated hard subprocess hh — )• H^W^ . 
III. ASSOCIATED H^W^ PRODUCTION 



Associated H^W^ production has two contributing subprocesses at leading order, quark- 
antiquark (gg) and gluon-gluon [gg) scattering. The leading order contributions correspond 
to the tree-level for gg, and one loop for gg fusion. Representative diagrams for these 
subprocesses are depicted in Figs. [1] and |2l respectively. Working in a five-flavor scheme 
(5FS) with an effective parton distribution for the h quark, the gg process is completely 
dominated by the hh contribution. Although the cross section for gg fusion is formally 
suppressed by two powers of the QCD coupling relative to hh annihilation, it may yield a 
comparable contribution at LHC energies due to the large gluon density at small x and needs 
to be taken into account. In the 5FS, there are additional contributions at higher orders of 
as where one gluon splits into a bb pair, giving hg — >■ H^W^h. These contributions should 
in principle be matched to the hh — )■ H^W^ process, which would yield an effective QCD 



correction factor slightly less than one [ij]. For simplicity we ignore this factor throughout 
this work. If the charged Higgs boson is light enough (m//± < rrit), there is an additional 
contribution to H^W^ production through top quark decays t — >■ hH^. When s is close to 
the tt threshold s ~ 4m^, on-shell top quarks can therefore give an additional production 
channel gg —> tt H^W^hb, which is enhanced by top resonances. However, in this 
work we study the case when the boson is heavy enough to be above the threshold for 
production in top decays, and this extra production channel is not relevant. 

The next-to-leading order (NLO) corrections to the H^W^ cross section in the MSSM are 
known (l9l-23|. but for the present study the leading order contributions suffice since we are 
mainly interested in comparing the NMSSM to the MSSM. We will, however, account for the 
most important higher order contributions from running quark masses, loop corrections to 
Higgs masses and mixings, and including (mass-dependent) widths of Higgs bosons appearing 
in s-channel propagators. The treatment of these effects is described in further detail below. 

Consider first the hh — )■ H^W^ contribution, which to the leading order is given by the 
tree-level diagrams in Fig. [T] At the parton-level, this contribution typically dominates over 
the gluonic one considered below. The corresponding parton-level cross section has the same 
form as in the MSSM case 
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m'* cot^ /3 ( 2m^ + P±) + fTT-t ^an"^ /3 ( 2m^p^ + t 
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where Gp is the Fermi constant, s, i, u are the usual Mandelstam variables, p±_ is the 
transverse momentum of the W boson in the hb cm. frame, and A(x, z) = + + — 
2{xy + yz + zx) is the Kallen function. The first line in Eq. fl3.1l) represents the s-channel 
resonance contribution, the last line corresponds to the non-resonant top quark exchange in 
the t-channel, and the second line contains the interference term. 

The functions Sq and Vq contain the propagators and relative couplings for the neutral 
Higgs bosons to quark flavor q. In Eq. (13.11) only the b quark contribution is needed, but for 
the gg contribution discussed below we need also the corresponding expressions for the top 
quark. These functions are defined as 



Si2 {Si2 COS/3 - 5'asin/3) 
sin/3 s — m?rT + imw- T 

1 S'ii (5*12 cos /3 - S'ii sin /3) 



COS 13 s-mjj^+ irriH, Th, 



1 ^ P.2(i"i2Cos/3-Pasin/3) 



sin /3 ^ ~ + r^, 



1 Pji {Pi2 cos /3 - Pg sin /3) 



where Th^ and F^. are the total (mass-dependent) decay widths of the and Ai bosons, 
respectively. We have obtained the expressions given in Eq. (13. 2p by modifying the S, V 
functions given in with the appropriate Yukawa couplings for the NMSSM case. We ne- 
glect the Yukawa couplings of the first- and second-generation quarks, as their contributions 
to the amplitude are negligibly small. If the masses of two (or more) neutral Higgs bosons 
with the same CP properties become degenerate, then the approximation used in Eq. f|3.2D 
breaks down, and one has to take into account Higgs mixing effects (see e.g. Refs. |48|, |49| ) . 
For the NMSSM scenarios we shall consider below, the masses will however always be such 
that Eq. (13. 2p remain valid. 

Let us now turn to the gg contribution. In analogy to the MSSM case 

the 

resonant amplitude of the gg — ?■ Hi^2,3, ^1,2 — ^ H^W^ subprocess from quark loops is given 
by the sum of all triangle diagrams of the type shown as the first diagram (upper left) in 
Fig. [2j This contribution can be written as 



V2 



Vxw = — as{fiR)GFmwe*Jpw){qi + ^2)^ 

71 



{dki - '^9"')m + ie^^'"^gipg2.H(S) 



(3.3) 



where ^^(/iR) is the strong coupling evaluated at the renormalization scale /xr, e* is the 
polarization vector of the W boson with momentum pw and helicity Xw, and are the 
polarization vectors of the gluons with momenta qi^2- These are summed over the color 
index c. The functions S and H come from the loop integration and correspond to neutral 
CP-even -^1,2,3 exchanges (S) and neutral CP-odd Ai^2 exchanges (H) in the s-channel. They 
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FIG. 2. Example leading-order diagrams contributing to the gluon-initiated hard subprocess gg — >■ 
H^W^. Dashed lines with arrows represent squarks, while internal dashed lines without arrows 
represent neutral Higgs bosons (-^1,2,3; ^1,2)- Only the diagrams with internal s-channel Higgs 
propagators will give rise to the resonant enhancement we are discussing in this paper. 



are given by 



s + ie 



(3.4) 
(3.5) 



where the sums run over all quark flavors q in the triangle loops, and the functions 



S{r) 



r 



1 — ( 1 1 arcsinh^V— r 

r 



P(r) = — arcsinh^ -y/^^, 
r 



(3.6) 
(3.7) 
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must be continued analytically for three regions in r, such that for r < 0,0 < r < 1, or 
r > 1, arcsinhA/— r must be represented by arcsinhA/— r, — iarcsiuA/r, or arccoshi/r — i7r/2, 
respectively. The contribution to the parton-level cross section is then given by 



do^ 

di 



■{gg ^ H+W~) 



20487r3 



A(s,m^,m^±)(|S(s)|2 + |n( 



(3.8) 



Due to Bose symmetry, the gg — >■ H^W^ cross section is symmetric with respect to t u 
interchange. Additionally, since we only consider the CP-invariant case, the cross sections 
for the gg — )■ H^W~ and gg — >■ H~W^ channels coincide. 

In our numerical calculations we take all possible quark and squark loop contributions 
into account from both triangle and box diagrams. For simplicity, we show here only the 
formulas for quark triangles, and do not list the complicated expressions for either the boxes 
or the diagrams with squark loops, schematically shown in Fig. |2l The full result also 
includes interference between these different contributions. We have checked our numerical 
results in the MSSM limit (which will be described below) against previous results from the 
literature 

Formally, the leading order contributions (as given by Egs. (|3.2p -f l3.8p above) contain 
tree- level masses and couplings. As advocated previously (see |2J] and references therein), 
higher order QCD and electroweak corrections can significantly affect MSSM observables. 
In particular the bottom Yukawa coupling is subject to large quantum corrections in the 
MSSM — as well as in the NMSSM — and these need to be taken into account properly. For 
this purpose, we follow the general recipe given in j2^. To take into account the large 



(SM) QCD corrections to the leading-order result, we use the QCD running b quark mass 



m 



DR 



(/xr). At two- loop order it is given by [50| 



37r 



a: 



1447r2 



(73 - 3n) 



(3.9) 



where n is the number of active quark flavors and (jj^ji) is the standard MS running mass 
4.2 GeV as input). Then, including the tan /3-enhanced supersymmetric 



we use m 



MS, 



nib 



QCD (SQCD) and electroweak (SEW) corrections |5ll] by a straightforward generalization 



of the MSSM results |52|, we obtain the following effective bottom-Higgs couplings: 



A 



off 
bbHi 



m 



DR 



Sil 



V2vcosf3l + Ab 



Ab 



5. 



i2 



X 



eff 
bbAk 



m 



DR 



V2vcosf3l + Ab 



1 + Ab 



Sil tan f3 

P^2 



1,2,3 



Pii tan f3 



k = l,2 



(3.10) 
(3.11) 



In a similar manner we also include the relevant corrections to the H^th vertex [53[. These 
so-called A;, corrections consist of two dominating parts, = A^*^'"° + A^^^, absorbing the 
leading SQCD and SEW corrections. In our case, the latter is dominated by the Higgsino- 
stop contribution A^^^ ^ A^^. In complete analogy to the MSSM case we therefore have 
Ab - 



^SQCD ^ ^j^g^g 



A SQCD 
^6 



2«s(<5) . nrf 2 2 

m-g /ieff tan [5 1{m~^^ , m^^ 



Q = i^h +^b2 +"^9)/3, 



Af 



mt 



IGtt^w^ sin (3 



At /ieff tan (3 1{m 



2 2 I l2^ 
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and finally 



/(a, b, c) = —-. — ( ab\n - + 6c In - + ca In - 

[a — b)[b — c){c — a) \ b c a 



Here rrig denotes the gluino mass, and m^^, mj. (z = 1,2) the sbottom and stop masses. 
From these expressions it is clear that the corrections could become large for either large 
values of fies and/or large tan/3. 

We shall refer below to the calculation at leading order, including the improvements 
discussed here (higher order corrections to Higgs masses and mixing, the Higgs widths in 
the s-channel Higgs propagators, the running m^, and the SUSY corrections to the bottom 
Yukawa couplings), as the improved Born approximation. 

To calculate the cross sections and perform the numerical computations, we have modified 
and extended the MSSM model file jli] of FeynArts [H] to contain the relevant NMSSM 
couplings and the necessary steps to use the improved Born approximation as discussed 



above. The parton-level amplitudes have been computed with FormCalc [56[ and inte- 
grated numerically. For the evaluation of the scalar master tree- and four-point integrals in 
the gluon contribution we have used the LoopTools library [s^. The Higgs mass spectra. 



mixing, couplings and decay widths have been calculated using NMSSMToOLS [4^,145 . 



IV. RESULTS 

A. NMSSM parameter dependence and benchmark scenarios 

As a first step, we investigate the impact on the parton-level cross sections ai,i^fj±-^y^ and 
d'gg^H^W'f of varying the NMSSM parameters. This information will be useful for defining 
the benchmark scenarios we are going to study in more detail below. We start from a generic 
NMSSM scenario with parameter values chosen as follows (this is what we will below refer 
to as Scenario A): 

A = 0.25, K = 0.25, Ax = -235 GeV, 
y4« = -150GeV, tan/3 = 10, 

and perform variations around these values. For each parameter point the partonic cross 
sections are evaluated as a function of a/J, using the improved Born approximation as 
described in Section UTTl 

To be able to study genuine NMSSM effects on the H^W^ process, we first want to 
compare our results to those obtained in the MSSM limit. We therefore look at the behavior 
of the partonic cross sections during the gradual transition from the NMSSM point defined 
by the parameter set given by Eq. (14.11) to the corresponding MSSM limitQ The results are 
shown in Fig.[3l A striking difference between the NMSSM and the MSSM is the presence of 
resonant enhancement of the partonic cross sections in the NMSSM. These resonances can 
be attributed to the heavy neutral Higgs poles s = mj^^ and s = , which for the default 
parameters have masses = 462.6 GeV and = 349.3 GeV (see Table [T] below) . Since 
in the MSSM it is generally true that mH,A < f^H^, these resonant contributions vanish 
in the MSSM limit. They are therefore an inherent feature of the NMSSM. As we will 



^ Wc remind the reader that the MSSM hmit is defined by taking A, k 0, while keeping the ratio k/X 
and all other parameters fixed to their respective values. 
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V?(GeV) 

FIG. 3. Parton- level cross sections for H^W^ production by bb (black) and gg (red) initial states for 
\=K = 0.25 (dotted), 0.1 (short-dashed), 0.05 (long-dashed), and in the MSSM limit X = k = W'^^ 
(solid). The remaining NMSSM parameters are fixed according to the description given in the text. 

show below, the resonant contributions are sensitive to the NMSSM parameters, and can 
give important contributions to the rate for pp — )■ H^W^ production. They are therefore 
interesting to study as a means of discriminating between the two models. 

Having established the presence of resonances as a potentially important difference for 
H^W^ production between the MSSM and NMSSM, we now proceed to study how the 
characteristics of these resonances are affected by variation of the NMSSM parameters. The 
results are shown in Fig. HI Starting from the upper left plot, we first consider different 
values of Ax. This has an immediate effect on the threshold through the linear dependence 
of rn'jj± on Ax (see Eq. (12. 3p ). It can also be seen that a change in Ax affects the presence 
of a resonance peak in the cross section. Qualitatively this can be understood from the fact 
that a charged Higgs boson with lower mH± can be produced from a (nearly) singlet scalar 
(pseudoscalar) of fixed mass, while this possibility disappears for a particular mass as mH± 
is increased. Further details on this point will be given below. Moving to the next plot, we 
see that a variation of A^ affects the cross sections in quite a different manner; it affects the 
relative positions and sizes of the resonance peaks. Effectively, when going from low \Ai^\ 
to higher values the two peaks, which at first are well separated, first meet and then get 
separated again, having changed positions. Since it also affects the relative size of the peaks, 
this signals that the mixture of the neutral Higgs resonances changes with modified A^. The 
two plots in the second row show the variation with A and k, respectively. In addition to 
the similar shift in threshold as demonstrated for the Ax case, we also see that these cross 
sections appear to depend on these parameters in the same combination as appears in mH± 
for fixed fies, thai is through the ratio n/X. On the last row of Fig. IHwe show the variation 
of the cross sections with /ies. This gives a very similar effect as choosing different values 
for Ax, which motivates us to consider a fixed value fies = 250 GeV in the following and 
instead use variation of Ax to control the value of mH± ■ The final plot shows the effect of a 
tan /3 variation. Also this variable enters the determination of (and thereby shifts the 
threshold), although in a more indirect way than fi^g or Ax- As can be seen from the figure, 
changing tan /3 also has a drastic effect on the absolute normalization of the cross sections 
and the width of the resonances. This comes mainly from the coupling of the charged Higgs 
boson to fermions of the third generation. For equivalent kinematic configurations, we can 
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FIG. 4. Parton-level cross sections for H^W^ production by bb (black) and gg (red) initial states 
versus \/I for different values of the NMSSM parameters, as indicated in the plots. The remaining 
NMSSM parameters are fixed according to the description in the text. 



therefore expect enhancements of the cross sections either for small or large values of tan /3. 

Based on these parameter variations we define five benchmark scenarios for further study 
(called scenarios A-E). The scenarios are selected to capture the different features of the 
partonic cross sections discussed above, and the values for the input parameters in the five 
benchmarks are shown in Table [H This table also gives the resulting Higgs mass spectrum, 
and the 5*3 3, P2,3 elements of the Higgs mixing matrices which give the singlet fractions of 
the heavy Higgs bosons H3 and A2. They will be relevant for the discussion of the resonant 
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Parameter 


Scenario 
ABODE 


Ax (GeV) 
(GeV) 
A 

K 

tan /3 


—235 —235 —235 —185 —243 
-150 -250 -400 -150 -150 
0.25 0.25 0.25 0.5 0.25 
0.25 0.25 0.25 0.5 0.25 
10 10 10 2.2 40 


Higgs mass spectrum (GeV) 




174.3 174.3 174.3 195.3 171.7 




118.4 117.4 115.0 114.6 120.3 

173.5 174.1 174.3 203.6 246.0 

462.6 435.3 391.1 459.6 463.3 


mAj_ 
ruA^ 


139.0 156.4 165.4 92.0 213.2 
349.3 438.2 549.2 383.4 355.7 


Singlet elements of H3, A2 


83,3 

P2,3 


0.993 0.991 0.986 0.988 0.992 
0.945 0.981 0.993 0.875 0.897 



TABLE I. Selected NMSSM benchmark scenarios, the corresponding Higgs mass spectrum, and 
singlet elements ^3^3, ^2,3 of the Higgs mixing matrices for the neutral heavy Higgs bosons. 



contributions below. 

Table [T] shows some features which are common to all the scenarios. As a general strategy 
we choose the parameters to obtain a rather low mH± in all scenarios (but still keeping 
mH± > rrit). For compatibility with LEP constraints 57|], we make sure that the light- 
est CP-even Higgs mass > 114.4 GeVH The scenarios are also compatible with the 
limits from direct Higgs searches at the Tevatron and the LHC implemented in HlGGS- 
BoUNDS v.2.1.0 [H^. Going into the specifics of the individual scenarios, the only difference 
between the definitions of scenarios A, B, and C is the value for A^. This variation leads 
to three different hierarchies for rriH^ and m^j, as can be read off from the table. Scenarios 
D and E are both mostly similar to Scenario A in terms of the heavy Higgs mass structure. 
Here we instead consider two "extreme" cases of low (D) and high (E) values for tan/5. 



B. Parton-level cross sections 

We now proceed to study associated H^W^ production in the NMSSM, making use of 
the benchmark scenarios defined in the previous section. The parton-level cross sections are 
evaluated as described in Section IIIIl and the results for ai,i^fj±^/^ are shown in Fig. O 

In this figure the solid lines give the NMSSM cross sections, while the dashed lines are 
the corresponding cross sections in the MSSM limit. As already discussed above, the most 

^ Applying this limit to the NMSSM in general is a very conservative approach, since it strictly speaking 
only applies to a Higgs with SM-like couplings. Specific NMSSM scenarios admit rriHi ^ 114.4 GcV 
without being in confiict with experimental data. However, in this work the selected scenarios correspond 
to the case where Hi is SM-like, and therefore the SM limit applies. 
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FIG. 5. Partonic cross sections o'^^ )H^w^ NMSSM benchmark scenarios A— E (solid lines), and 
the corresponding results in the MSSM limit (dashed). 



striking difference is the presence of resonances in the NMSSM case. In all the five bench- 
marks the resonances are manifest as either one or two peaks of varying size. Comparing the 
plots for scenarios A-C, we observe a shift of the large peak from low energies (for Scenario 
A) to high energies (Scenario C). There is also a corresponding down-shift of the smaller 
peak to low energies. For Scenario B the two peaks nearly coincide. The positions of these 
poles are determined by rriH^ and , and we note that the largest peak corresponds to the 
A2 resonance. This means that even if the singlet component of A2 is larger than that of H^, 
(1-^2,31 > I'S's.sl), the A2 couples more strongly to the hh initial state. For Scenario D we only 
observe one (small) resonance peak. This is due to the low value of tan/3, corresponding to 
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FIG. 6. Partonic cross sections o'gg_j>/f±vi/T for NMSSM benchmark scenarios A-E (solid lines), 
and the corresponding results in the MSSM limit (dashed). 



a reduced A2hh coupling. On the other hand, we also observe a larger overall cross section 
in the continuum (visible also in the MSSM limit), since the non-resonant contribution me- 
diated by t-channel top exchange increases as tan /3 <^ 7. For Scenario E (which has a high 
value of tan (5) something similar is observed for the continuum, but with an even larger 
enhancement of the cross section at low energies and a more rapid drop as — oo. In this 
case we also observe two pronounced resonances (as in scenarios A and C), but the fairly 
large couplings in Scenario E lead to larger differences between the NMSSM results and the 
MSSM limit also for energies away from the actual poles. 
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Turning now to (Tgg^H±w^^ the contribution from gluon fusion is expected to be richer 
than that initiated by 6-quarks, since it involves additional non-resonant box diagrams. 
The interference with these can strongly affect the resulting cross section (for a study of 
these interference effects in the MSSM, see The cross section is again evaluated as 

outlined in Section IIIII and the results are shown in Fig. El One look at this figure reveals 
that the gluon-initiated process has a much lower cross section compared to the hb initial 
state, about 3-4 orders of magnitude. We can however expect this difference to be (at 
least partly) compensated in the hadronic cross section by the larger gluon content of the 
proton at intermediate and small x (see below). Compared to the hh process, we do observe 
a general broadening of the resonances, and larger differences between the NMSSM and 
the MSSM limit — the latter in particular for energy ranges between near-lying resonances, 
where interference can lead to either an enhanced or suppressed cross section prediction in 
the NMSSM compared to the MSSM. Most of the gg distributions show a feature at the 
top pair threshold = 2mt (sometimes masked by an NMSSM resonance). Since this 
kinematic effect is present also in the MSSM it is not interesting for the comparison of the 
results between the two models. 

Looking specifically at the results for scenarios A-C in Fig. |6l they display the same 
resonance structure as the hb case. This tells us that both the and A2 resonances 
play a role also here. However the peaks are more similar in size (for Scenario A), and in 
the case of Scenario C we see that the low energy peak is dominating the cross section. 
Since this corresponds to the contribution, we conclude that the resonant process here 
is instead dominated by the H^ti coupling which enters the top loop contribution. The 
importance of the coupling to the top for lower tan/3 becomes evident for Scenario D. This 
scenario has a greatly enhanced cross section, both compared to scenarios A-C, and with 
very pronounced resonance enhancements compared to the MSSM limit. Also Scenario E 
benefits from the same large continuum cross section as observed for the hh contribution, 
but with the additionally boosted resonance contribution of observed generally for the 
gg process. 



Hadron- level cross sections 



The total hadronic cross section (Jpp^H±y/T is obtained from the partonic cross sections 
dij (with ij = hb,gg) by integration over the parton distribution functions (PDFs). This 
can alternatively be expressed in terms of parton luminosities, which allows studying the 
impact of the PDFs on the cross section. The parton luminosities Lij for partons i,j are 
defined as 



dLij 
dr 



where t = X1X2 and fi{x, fip) is the PDF for parton i evaluated at the factorization scale /ip. 
Note that, since we are interested in bb and gg scattering, we will always have fi{x) = fj{x). 
The significance of r is that the center-of-mass energy of the partonic system is given by 
= i/ri, where ^/s = Eqm is the collider energy. Using the parton luminosities, the total 
cross section can be expressed as 



'pp- 



dr 


'1 dLij' 


/o ^ 


s dr 



[saijis)] 



(4.3) 
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Note that the factor involving the parton luminosity has the dimensions of a cross section; it 
can be used to estimate the size of the hadronic cross section when the partonic cross section 
is known. Even though the gg cross sections are much smaller than the hb cross sections at 
parton-level, at hadron-level the much larger gluon PDFs make the gg contribution compet- 
itive with the hh contribution. In Fig. [7] we show the ratio of the gg and hh luminosities at 7 



TeV and at 14 TeV, calculated with MSTW PDFs [59| . We see see from this figure that the 
gg/bb luminosity ratio is typically a factor 1000 at 7 TeV, and slightly smaller for 14 TeV. 

For the calculation of the hadronic cross sections, we modify the benchmark scenarios 
listed in Table |T] to allow for variation of mH±- This is achieved by varying the value of Ax, 
keeping the other parameters fixed to the values given in the table. Besides changing mH±, 
this affects the doublet-dominated neutral Higgs bosons (which never appear resonantly), 
while the masses of the singlet-dominated Higgses {H^ and A2 for low mH±) are largely 
insensitive to Ax variations. This means in particular that the resonance structure discussed 




FIG. 8. Dependence of the neutral Higgs masses on mf^± in Scenario A with varying charged Higgs 
boson mass. The different curves correspond to (in black): tjihi (sparse dots), mj^j (dotted), niHs 
(dense dots), and in green rriA^ (dot-dashed, short), and m^^ (dot-dashed, long). The solid red 
line shows the threshold for associated H^W^ production. 
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above will remain unaltered over certain ranges for m//± (when the doublet mass is well below 
the singlet mass scale). This important point is demonstrated in Fig. |8l which shows the 
dependence of the five neutral Higgs masses in Scenario A on mH± (when varying Ax). A 
qualitatively similar picture is obtained in the other scenarios we study. 

For the numerical evaluation of the cross sections we consider pp collisions at the LHC 
at the two center of mass energies ^/s = 7 TeV and a/s = 14 TeV. We use CTEQ6.6 
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FIG. 9. Total hadronic cross section pp — ?• H'^W^ at = 7 TeV (left) and ^/s = U TeV (right) 
as a function of the charged Higgs boson mass in scenarios A~C (with variable ^a)- The symbols 
indicate the separate contributions to the NMSSM cross section from bb (dots) and gg (crosses), 
while the lines show the corresponding contributions {bb solid, gg dotted) in the MSSM limit. 
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FIG. 10. Total hadronic cross section pp — ?• H^W^ at ^/s = 7 TeV (left) and = U TeV (right) 
for NMSSM Scenario D. The different symbols used are explained in the caption of Figure [H 



parton distributions [OOj and a fixed fip = + m\y- The results for Scenarios A-C are 
presented in Fig. [9l shown in parallel for 7 TeV (left) and 14 TeV (right). We give the 
contributions of bb (big dots) and gg (crosses) separately, together with their respective 
contributions in the MSSM limit (as solid and dotted lines, respectively). By comparing 
the NMSSM results to the MSSM limit, it can be seen from these plots that it is possible 
to have substantial NMSSM enhancements also in the hadron-level cross section. In some 
mass ranges this enhancement can be an order of magnitude. The observed shapes, with an 
almost fiat dependence of the cross section on the charged Higgs mass up to some mass, and 
then a rapid fall-off to the MSSM value, is caused by the contribution of the resonances. In 
the region where the H^W^ threshold is below the resonance mass (cf. Fig. |8]), the cross 
section is dominated by the contribution of resonant diagrams. When the transition takes 
place, the non-resonant behavior quickly becomes similar to the MSSM case. Note that the 
contributions from bb and gg in general receive their largest resonant enhancements from 
two different resonances, as discussed in Section IIVBI This also leads to the two different 
mfj± regions where to cross section is flat for the two cases. Similarly to the discussion of 
the partonic cross sections above, we can see how the three scenarios A-C differ mainly in 
the position of the resonances, which results in a cross section enhancement (for bb) at low 
mH± in Scenario A, and towards higher masses in Scenario C. The opposite is true for the 
gg contribution. 

In Fig. [TU]we show the hadronic cross section for Scenario D. We see from this figure that 
the low value of tan /3 in this scenario gives a significantly larger gg contribution to the total 
cross section than in the other cases; it dominates over bb for the full range oi mH±. There 
is more room for a large gg contribution in the NMSSM, since lower tan /3 is still allowed in 
this model compared to the MSSM. We note a relatively small resonance enhancement in 
this scenario, in particular for bb where it is barely visible. For our last scenario. Scenario 
E, the cross section is shown in Fig. [TTl With a large tan /3, this scenario has been selected 
to maximize the Higgs couplings to bb, and indeed we find the largest total cross section for 
this scenario. As could be expected, the bb contribution dominates completely, being almost 
two orders of magnitude larger than gg for mj^± > 400 GeV. In this scenario the resonances 
are again pronounced, and we observe clearly the coupling of bb to the (lighter) A2, while 
the gg contribution is more efficiently enhanced by the (heavier) resonance. 

It is interesting to note from the hadron-level results that the relative importance of bb 
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FIG. 11. Total hadronic cross section pp — t- H^W^ at ^/s = 7 TeV (left) and = U TeV (right) 
for NMSSM Scenario E. The different symbols used are explained in the caption of Figure [9j 
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FIG. 12. The gluonic contribution to the total hadronic cross section pp — )• H^W^ at ^/s = 7 TeV 
(left) and ^/s = U TeV (right) for NMSSM Scenario A (crosses) compared to the MSSM hmit 
(dotted line). The blue lines correspond to the separate contributions to gg — )■ H^W^ from only 
including triangles (solid line) and box diagrams (dashed). 



and gg may change drastically in regions where either of the resonances dominate. This 
change in the relative contributions can be much larger than what is expected from the 
energy scaling of the parton luminosities, cf. Fig. [71 It is a result of the different couplings 
of the and A2 resonances to the initial and final states. In some of our scenarios the gg 
contribution can even become larger than the bb contribution in certain regions of mH± (in 
Scenario D, the neutral Higgses have a reduced coupling to b quarks, and the gg contribution 
is dominant for all values of mH±). 

Another feature of the gg contribution to the hadronic cross section which deserves a 
comment is the apparent larger NMSSM enhancement at the 7 TeV LHC compared to the 
14 TeV case. The shape of the enhancements is due to a convolution of the resonance 
peaks observed in the parton-level cross section with the dependence of the PDFs on the 
partonic center-of-mass energy \f§ (including the xi^2-dependence of the PDFs). There can 
also be interference effects between the resonant contributions and (non- resonant) boxes. 
To investigate the relative importance of the different contributions in some more detail, we 
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illustrate in Fig. [T2]tlie contributions from (boxes)^ and (triangles)^ to the total cross section 
for Scenario A. At y/s = 7 TeV, we see that the box and triangle contributions are similar 
in the non-resonant region, i.e. at relatively large charged Higgs masses mj^± > 380 GeV. 
Interestingly enough, the full cross section is numerically at the same level as separately the 
box and triangle contributions, which means that the destructive interference between the 
two is large and similar to that in MSSM. However, in the resonant region (corresponding to 
< 380 GeV), the triangle contribution becomes much more pronounced and strongly 
enhances the total cross section. The interference effects are naturally quite small there. At 
higher energies, ^/s = 14 TeV, we observe a somewhat different picture in Fig. [121 The boxes 
in this case gives a large (dominant) contribution to the total gluon-initiated cross section 
over the whole range inmH±- In the non- resonant region, m//± > 380 GeV, the interference 
with the smaller triangle contribution noticeably decreases the cross section compared to the 
box contribution alone. In the resonant region, mH± < 380 GeV, the triangles become more 
important, but remain sub-dominant compared to the boxes. The cross section is therefore 
only enhanced slightly with respect to the MSSM case. Analogously to the 7 TeV case, the 
interference turns out to be less important than in the non-resonant region. 



V. COLLIDER PHENOMENOLOGY 

We have shown that the total cross sections for associated H'^W^ production in the 
NMSSM can be substantially enhanced compared to the corresponding cross sections in the 
MSSM. We would now like to discuss the phenomenological implications of this for searches 
at the LHC. Dedicated collider studies of this channel have been performed for MSSM and 



two-Higgs doublet models [22], |25|-|28|. While we will comment here on what features of 
these studies that are relevant for NMSSM — and what features will be different — we leave 
a dedicated phenomenological study for the future, since this would require calculation of 
differential cross sections and consideration of backgrounds and sensitivities. 

The defining feature of the various possible ways of detecting associated production is 
how the charged Higgs decays. The previous studies have all concentrated on the decays 
that are relevant in the MSSM or 2HDM, namely t^u,cs for light charged Higgs 

bosons and — tb for heavy charged Higgs. Decays to SUSY particles may also be 
important for heavy H^. As we discussed above, in the NMSSM the decay — )■ W^Ai 
is sometimes dominant, and will lead to quite different experimental requirements. This 
channel is normally not possible in the MSSM, since and mfj± are close to degenerate. 
In the NMSSM, the decay width is proportional to the doublet component of Ai, but may 
be large even if the Ai is mostly singlet [lo[ lH- The Ai boson can be very light in the 



NMSSM, and in such a case its dominant decays will be Ai — bb or Ai — )■ r+r". We do not 
consider scenarios with a light Ai in this paper, but they could nevertheless be of interest. 

In Fig. [T3]we show the decay branching ratios of H~^, calculated using NMSSMToOLS, 
for three of our benchmark scenarios. Since the most important parameter entering the 
determination of the decay modes is tan/3, the results for scenarios B and C are very similar 
to those for Scenario A and we do not show them explicitly. All our scenarios have, as 
explained above, charged Higgs masses above the threshold for the decay — > tb. As can 
be seen from Fig. [T31 the decay — )■ tb dominates over a wide range in masses, but for 
scenarios A and E the decay t^u is appreciable over the entire mass ranges plotted, 

and dominates close to the tb threshold. The — )■ ly^Ai decay is only relevant for 
Scenario D, which has a somewhat lighter Ai, or for very heavy in scenarios A-C. SUSY 
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FIG. 13. Charged Higgs branching ratios — )• t'^u (sohd, red), — )■ tb (long dashes, black), 
H+ cb (dotted, purple), H+ W'^Ai (short dashes, blue), and i?+ SUSY (dot-dashed, 
gray) for NMSSM benchmark scenarios A (left), D (center), and E (right). 



decays to a chargino-neutralino pair, — )■ xtXi^ also become appreciable above threshold, 
but this of course depends very much on the parameters entering the neutralino and chargino 
sectors. To summarize, the situation for our scenarios is thus similar to the MSSM case in 
terms of decay channels, except possibly for Scenario D. However, the possible Higgs masses 
and the size of the cross sections can be different. We also would like to emphasize that, for 
a complete coverage of the phenomenological possibilities in the NMSSM, other scenarios 
where the decay modes are different should be considered. 

Previous collider studies, which were all based on the MSSM, have thus used either the 
— )■ T^v or the — )■ th decay channels. We will now make some simple estimates of 
how these results translate to the NMSSM. Keep in mind that although we can draw some 
general conclusions based on these estimates, there may be more subtle differences, caused 
e.g. by different dependences of the differential cross sections on kinematical variables when 
the relative contributions of different types of processes are different. 

In [26|, the — >■ t^u channel was considered with hadronic r and decays. The 
main background to this signal is the production of 1^ + 2 jets. A parton-level study showed 
that for m//± = 175 GeV in the m™^^ scenario at tan /3 = 50, it should be possible after 
a series of cuts to get a statistically significant signal at the 14 TeV LHC with 300 fb~^ of 
data. In this scenario, the total production cross section was 55 fb, with a charged Higgs 
branching ratio into t^v of 100%, and the obtained S/\fB value was 17. We may rescale 
this result to e.g. our Scenario A (with a much lower tan/3), where the cross section is, 
conservatively, 100 fb and the branching ratio BR(if+ — t- t^v) = 100%, which would yield 
a significance of about S/^/B = 30. A scenario with m/^± = 400 GeV was also studied (26j . 
but here the significance was only S/ a/B = 2.6 due to the lower cross section and branching 
ratio. At this higher mass, we do not predict an enhancement of the cross section in the 
NMSSM for Scenario A, but rather a similar cross section to the one in the MSSM; our 
expected significance should therefore be similar. If we instead consider Scenario C, which 
has the resonant enhancement at a higher mass, the cross section at mH± = 400 GeV is 
10 fb, while the branching ratio is about 10% (instead of 20%). This would yield a higher 
significance, but still not above S/ \/B = 5. Similar conclusions as those presented in 26 
were found inj22[, and more recently in [23], where a detailed study of discovery contours in 
the MSSM was performed using both leptonic and hadronic W decays. 

The — )■ tb decay was first studied in [25| , where hadronic decays of the top quark and 
the W from the top quark, and leptonic decay of the other W were considered. The hadronic 
decay of the has the advantage over the tu channel that the mass can be completely 
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reconstructed, while the charged lepton from the other W provides a useful trigger. The 
conclusion of this work was unfortunately that the tt background completely overwhelms 



the signal. However, it was recently argued [28| that with the use of optimized cuts this 



channel can be useful to obtain significant results. No numbers for the significance of the 
— !■ tb channel have been given in j25|, but we note that in our Scenario A, the total 
cross section times branching ratio at mH± = 200 GeV is roughly 10 times larger than the 
one used in the analysis there, and it is therefore plausible that the significance could be 
improved. 

The remaining, possibly useful, decay channel is — )■ W^Ai, followed by Ai — > bb or 
Ai — r+r" decay. Note that the final state in the bb case is the same as for the hadronic 
— !■ tb case, but with a bb pair that should reconstruct the Ai mass. This may provide 
an additional handle on the signal, but it is not obvious that this is useful experimentally; it 
may be that the r+r" decay proves more useful. If possible, we would like to suggest to the 
experiments that tt samples with leptonic W decays are investigated for bb resonances. To 
look for (low mass) bb resonances in more energetic events has also been proposed recently as 
a means of searching for light, singlet-dominated, NMSSM Higgs bosons produced in SUSY 
cascades j6l|. In any case, we think that the — >■ W^Ai channel deserves a detailed 
study, and as far as we know this has not been performed for the case of H^W^ production. 



VI. SUMMARY AND CONCLUSIONS 

We have studied associated charged Higgs and W boson production in the NMSSM. 
This process is complementary to the main production modes anticipated for heavy charged 
Higgs bosons {mH± > frit) at the LHC We calculated the leading order contributions to the 
total hadronic cross section pp — )■ H^W^ in a general NMSSM setting, corresponding to 
the tree-level bb contribution and the g'g'-initiated subprocess at one loop. The calculation 
has been performed using an improved Born approximation, taking into account the most 
important effects of (S)QCD higher order corrections. 

For H^W^ production in the NMSSM, we have first investigated the parameter depen- 
dence of the parton-level cross section and the corresponding result in the MSSM limit. 
Starting from a maximal mixing scenario, we have then defined five NMSSM benchmark 
scenarios, with the common feature that they allow for resonant contributions from heavy 
singlet-like Higgs bosons. These resonances couple both to the bb and gg initial states, and 
we find that they can lead to significant enhancements of the cross section (by up to an 
order of magnitude) over a wide range of charged Higgs masses in the benchmark scenarios. 
The presence of these resonances is a genuine feature of the NMSSM, since Higgs mass con- 
figurations like this are not possible in the MSSM. This process might therefore be useful as 
a discriminator between the two models. 

We also discussed briefly the phenomenological implications of H^W^ production in the 
NMSSM. Based on previous work for the MSSM, we estimated the discovery significance 
that could be expected for different decay channels of H^. From these estimates it seems 
likely that the chances of detecting the charged Higgs boson of the NMSSM at the LHC 
through this process are quite good, at least in scenarios similar to ours with a not too heavy 
H^. It may not even be necessary to wait for the 300 fb~^ at 14 TeV that were previously 
assumed in MSSM studies, but clearly it will still require more data than "standard" Higgs 
searches. However, it should be remembered that these results are rough estimates based 
on rescaling of earlier MSSM results with the total cross section ratio. More detailed Monte 
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Carlo studies should therefore be done, for example on the impact of the different kinematics 
between the resonant and non-resonant contributions. It would also be necessary to revise 
previous results for the present and future running conditions of the LHC, something which is 
also true for the existing MSSM studies. To more seriously assess the prospects of observing 
associated H^W^ production at the LHC it would eventually be necessary to perform an 
analysis at the level of full experimental detector simulation. Perhaps the encouragement 
provided by these positive results could trigger the interest necessary to revisit this channel 
as a probe for the NMSSM and other models beyond the MSSM. 
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